	PRO MHDFIGRATE,filename1,filename2,filename3,filename4,figname
; 
;
loadct, 5
restore,filename1+'_rec_flux.idl'
t1 = tiempo
r1 = recflux
restore,filename2+'_rec_flux.idl'
t2 = tiempo
r2 = recflux/sqrt(0.001)
restore,filename3+'_rec_flux.idl'
t3 = tiempo
r3 = recflux/sqrt(0.005)
restore,filename4+'_rec_flux.idl'
t4 = tiempo
r4 = recflux/sqrt(0.01)

maxrec = max([max(r1),max(r2),max(r3),max(r4)])
maxtime = max([max(t1),max(t2),max(t3),max(t4)])

nu = '!4' + String("147B) + '!X'
set_plot,'ps'
device,filename=figname+'.eps',bits_per_pixel=8,/encapsulated,xs=15,ys=12,/color
plot,t1,r1,thick=3,charsize=1.5,xtitle='Time',tit='Reconnected Flux',xr=[0,10],yr=[0,maxrec],xsty=1,ysty=1
xyouts,0.5,0.27,nu+' = 0',charsize=1.5
oplot,t2,r2,thick=3,color=45
xyouts,0.5,0.25,nu+' = 0.001',color=45,charsize=1.5
oplot,t3,r3,thick=3,color=90
xyouts,0.5,0.23,nu+' = 0.005',color=90,charsize=1.5
oplot,t4,r4,thick=3,color=135
xyouts,0.5,0.21,nu+' = 0.01',color=135,charsize=1.5
device,/close
set_plot,'x'
end

